Int. J. Adv. Sci. Eng. Vol.9 No.2 2709-2735 (2022) 2709 E-ISSN: 2349 5359; P-ISSN: 2454-9967 


Optimization of Voltage-Controlled Oscillator based on 
Heterostructure to increase the density of their elements. 
Influence of Mismatch-Induced Stress and Porosity of 
Materials on Technological Process Development 


E.L. Pankratov 
Nizhny Novgorod State University, 23 Gagarin Avenue, 
Nizhny Novgorod, 603950, Russia 


ABSTRACT: In this paper we introduced an approach to increase integration rate of field-effect 
heterotransistors in the framework of a voltage controlled oscillator. We formulate recommendations 
for optimization of annealing to decrease dimensions of field-effect heterotransistors and to increase 
their density. At the same time one can decrease transition area between doped and undoped areas of 
elements of the considered oscillator. Also we consider an approach to decrease the value of 
mismatch-induced stress in the considered heterostructure. We introduce an analytical approach to 
model diffusion and ion types of doping with account concurrent changing of parameters in space and 
time. At the same time the approach gives us possibility to take into account nonlinearity of 
considered processes. 
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INTRODUCTION 

In the recent time several actual problems of 
the solid state electronics (such as increasing of 
performance, reliability and density of elements 
of integrated circuits: diodes, field-effect and 
bipolar transistors) are intensively solving [1- 
6]. To increase the performance of these devices 
itis attracted an interest for the determination 
of materials with higher values of charge 
carrier’s mobility [7-10]. One way to decrease 
the dimension of elements of integrated circuits 
is manufacturing them in thin film 
heterostructures [3-5,11]. In this case, it is 
possible to use inhomogeneity of 
heterostructure and necessary optimization of 
doping of electronic materials [12] and 
development of epitaxial technology to improve 
these materials (including analysis of mismatch 
induced stress) [13-15]. An _ alternative 
approaches to increase the dimensions of 
integrated circuits are using of laser and 
microwave types of annealing [16-18]. 


In the present study, we introduce theoretical 
approach will be useful for the manufacturing of 
field-effect transistors. The approach provides a 
possibility to decrease their dimensions with 
increasing their density, framework a voltage 
controlled oscillator. We also consider the 


possibility to decrease mismatch-induced stress 
to decrease quantity of defects, generated due 
to the stress. We consider a heterostructure, 
which consist of a substrate and an epitaxial 
layer (see Fig. 1), also consider a buffer layer 
between the substrate and the epitaxial layer. 
The epitaxial layer includes several sections, 
which were manufactured by using another 
material. These sections have been doped by 
diffusion or ion implantation to manufacture 
the required types of conductivity (p or n). 
These areas became sources, drains and gates 
(see Fig. 1). After doping, it is required to anneal 
of dopant and/or radiation defects. Main aim of 
the present paper is analysis of redistribution of 
dopant and radiation defects to determine 
conditions, which correspond to decreasing of 
elements of the considered oscillator and at the 
same time to increase their density. At the same 
time we consider a possibility to decrease 
mismatch-induced stress. 


METHOD OF SOLUTION 

To solve our aim we determine and analyzed 
spatio-temporal distribution of concentration of 
dopant in the considered heterostructure. We 
determine the distribution by solving the 
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second Fick's law in the following form [1,20,23]. 
Ley Ea 
cet i 


Fig. 1a Structure of the considered oscillator [19] 
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Substrate 


Fig. 1b Heterostructure with a substrate, epitaxial layers and buffer layer (view from side) 


eel nel) 2 | p2Clennt)), 2| poClinel)|, 2| pockn2)), 
Ot Ox Ox Oy Oy Oz Oz 


+ 08 Pev.u(s eee. yw) + 
Ox| kT 0 


L, 
+02 | Pov asils y,z,t) (C(x, pwa)aw | (1) 
Oy| kT 0 
+O] Pes 6 wley.zt)|, 0 | Des 2 wleyzt)), 2) Des 2 w(ey,z0) 
Ox|VkT Ox Oy|VkT Oy Oz|VkT OZ 


with boundary and initial conditions 


AC (x, y,z,t) 
Ox 


_0 AC (x, y, Zt) 
Ox 


éC(x, y, Zt) 


=0, 
Oy 


= 0 ,C (x,y,2,0)=fc (x,y,Z), 
y=0 


x=0 x=L, 


Pankratov., 


International Journal of Advanced Science and Engineering www.imahendrapublications.com 


Int. J. Adv. Sci. Eng. Vol.9 No.2 2709-2735 (2022) 2711 


oC (x, ae t) 

Oy Ls Oz 

Here C(xy,z,t) is the spatio-temporal 
distribution of concentration of dopant; © is the 
atomic volume of dopant; V; is the symbol of 


L, 
surficial gradient; [Cc (x, y,z,t)dz is the 
0 


surficial concentration of dopant on interface 
between layers of heterostructure (in this 
situation we assume, that Z-axis is 
perpendicular to interface between layers of 
heterostructure); s4(x,y,z,t) and se(x,y,z,t) are 
the chemical potential due to the presence of 


D, =D,(x,y,z,T)|1+é . 


C 
D; =D (a; y> z,T) 1+és 


Here D; (x,y,z,T) and Dis (x,,z,T) are the 
spatial (due to accounting all layers of 
heterostruicture) and temperature (due to 
Arrhenius law) dependences of dopant diffusion 
coefficients; T is the temperature of annealing; P 
(xy,z,T) is the limit of solubility of dopant; 
parameter y depends on properties of materials 
and could be integer in the following interval y 
€[1,3] [24]; V (x,y,z,t) is the spatio-temporal 
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mismatch-induced stress and_ porosity of 
material; D and Ds are the coefficients of 
volumetric and surficial diffusions. Values of 
dopant diffusions coefficients depends on 
properties of materials of heterostructure, 
speed of heating and cooling of materials during 
annealing and spatio-temporal distribution of 
concentration of dopant. Dependences of 
dopant diffusions coefficients on parameters 
could be approximated by the following 
relations [24-26] 
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distribution of concentration of radiation 
vacancies; V’ is the equilibrium distribution of 
vacancies. Concentrational dependence of 
dopant diffusion coefficient has been described 
in details in [24]. Spatio-temporal distributions 
of concentration of point radiation defects have 
been determined by solving the following 
system of equations [20-23,25,26] 
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Here I! (x,y,z,t) is the spatio-temporal 
distribution of concentration of radiation 
interstitials; /* is the equilibrium distribution of 
interstitials; D)(x,y,z,T), Dv(xy,z,T), Dis(xy, Z,T), 
Dys(x,y,Z,T) are the coefficients of volumetric 
and surficial diffusions of interstitials and 
vacancies, respectively; terms V(x,y,z,t) and 
[2(xy,z,t) correspond to generation of 
divacancies and diinterstitials, respectively (see, 
for example, [26] and appropriate references in 
this book); kiv(xy,z,T), kixy,z,T) and 
kvyv(x,y,Z,T) are the parameters of recombination 
of point radiation defects and generation of 


kT ia + y, +z, 


their complexes; k is the Boltzmann constant; @ 
= a3, a is the interatomic distance; ¢ is the 
specific surface energy. To account porosity of 
buffer layers we assume, that porous are 
approximately cylindrical with average values 


x =F y; and z; before annealing [23]. 


With time instead small pores one can find 
series vacancies. Modification of distribution of 
vacancies in time shows absorption of the above 
vacancies by the above pores [27]. With time 
large pores became larger due to absorbing the 
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vacancies and became more spherical [27]. heterostructure, existing due to porosity, could 
Distribution of concentration of vacancies in be determined by summing on all pores, ie. 
lm n . ; 7 ; 5 
Vixny2t=x>d > V,(x+ia,yt+jB.zt+k x,t), R=yx +y +z. 
i=0j=0k=0 


Here a, Zand yare the average distances between centers of pores in directions x, y and z; 1, mandn 
are the quantity of pores inappropriate directions. 


Spatio-temporal distributions of divacancies @y (x,,z,t) and diinterstitials ® (x,y,z, t) could be 
determined by solving the following system of equations [25,26] 
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Here Do(x,y,z,T), Dov(Xy,Z,T), Dor (Xy,Z,T) and Dovs(x,y,z,T) are the coefficients of volumetric and 
surficial diffusions of complexes of radiation defects; ki(x,y,z,T) and ky(x,y,z,T) are the parameters of 
decay of complexes of radiation defects. 


Chemical potential z in Eq.(1) could be determine by the following relation [20] 
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tensor; uj uj; are the components u;(x,y,z,t), uy(xy,z,t) and u,(x,y,z,t) of the displacement vector 
Uu (x, y,Z, t); Xj, Xj are the coordinate x, y, z. The Eq. (3) could be transform to the following form 
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Integration of the left and right sides of the Equation (1b), (3b) and (5b) on time gives us possibility to 
obtain relations for above approximation in the final form. 
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We determine average values of the first-order approximations of concentrations of dopant and 
radiation defects by the following standard relation [28] 
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We determine the approximations of the second n-th order of concentrations of dopant and 
and higher orders of concentrations of dopant radiation defects we replace the required 
and radiation defects framework standard concentrations in the Eqs. (1c), (3c), (5c) on the 
iterative procedure of method of averaging of following sum Qnpt pn-1(X,y,Z,t). The replacement 
function corrections [28]. Framework this leads to the following transformation of the 
procedure to determine approximations of the appropriate equations 
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We determined the concentration of dopant, 
concentrations of radiation defects and 
components of displacement vector by using the 
second-order approximation framework 
method of averaging of function corrections. 
This approximation is usually enough good 
approximation to make qualitative analysis and 
to obtain some quantitative results. All obtained 
results have been checked by comparison with 
results of numerical simulations. 


DISCUSSION 

We have analyzed the dynamics_ of 
redistributions of dopant and radiation defects 
during annealing and under influence of 
mismatch-induced stress and modification of 
porosity. Typical distributions of concentrations 
of dopant in heterostructures are presented on 
Figs. 2 and 3 for diffusion and ion types of 
doping, respectively. These distributions have 
been calculated for the case, when value of 
dopant diffusion coefficient in doped area is 
larger, than in nearest areas. The Figures show 

LyL, 


that the inhomogeneity of heterostructure gives 
us possibility to increase compactness of 
concentrations of dopants and at the same time 
to increase homogeneity of dopant distribution 
in doped part of epitaxial layer. However, this 
approach of manufacturing of bipolar transistor 
it is necessary to optimize annealing of dopant 
and/or radiation defects. Reason of this 
optimization is as follows: If annealing time is 
small, the dopant did not achieve any interfaces 
between materials of heterostructure. In this 
situation one cannot find any modifications of 
distribution of concentration of dopant. If 
annealing time is large, distribution of 
concentration of dopant is too homogenous. 
Recently, we optimized annealing time [29-37]. 
Framework this criterion we approximate real 
distribution of concentration of dopant by step- 
wise function (see Figs. 4 and 5). Farther we 
determine optimal values of annealing time by 
minimization of the following mean-squared 
error 


ic (x, y,z,0)-y (x,y, z)]d zd ydx, (15) 
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Fig.2. Distributions of concentration of infused dopant in heterostructure from Fig. 1 in direction, 
which is perpendicular to interface between epitaxial layer substrate. Increasing of number of curve 
corresponds to increasing of difference between values of dopant diffusion coefficient in layers of 
heterostructure under condition, when value of dopant diffusion coefficient in epitaxial layer is larger, 
than value of dopant diffusion coefficient in substrate 
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Fig.3. Distributions of concentration of implanted dopant in heterostructure from Fig. 1 in direction, 
which is perpendicular to interface between epitaxial layer substrate. Curves 1 and 3 corresponds to 
annealing time © = 0.0048(L,2+L,2+L,2)/Do. Curves 2 and 4 corresponds to annealing time © = 
0.0057(L,2+Ly2+L72)/Do. Curves 1 and 2 correspond to homogenous sample. Curves 3 and 4 
corresponds to heterostructure under condition, when value of dopant diffusion coefficient in 
epitaxial layer is larger, than value of dopant diffusion coefficient in substrate 
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Fig. 4. Spatial distributions of dopant in heterostructure after dopant infusion. Curve 1 is idealized 
distribution of dopant. Curves 2-4 are real distributions of dopant for different values of annealing 
time. Increasing of number of curve corresponds to increasing of annealing time 


Where y (x,y,z) is the approximation function. necessary to anneal radiation defects after ion 
Dependences of optimal values of annealing implantation. One could find spreading of 
time on parameters are presented on Figs. 6 and concentration of distribution of dopant during 
7 for diffusion and ion types of doping, this annealing. In the ideal case distribution of 
respectively. It should be noted, that it is dopant achieves appropriate interfaces between 
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materials of heterostructure during annealing of 
radiation defects. If dopant did not achieve any 
interfaces during annealing of radiation defects, 
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dopant. In this situation optimal value of 
additional annealing time of implanted dopant 
is smaller, than annealing time of infused 
dopant. 


it is practicably to additionally anneal the 
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Fig. 5. Spatial distributions of dopant in heterostructure after ion implantation: Curve 1 is idealized 
distribution of dopant. Curves 2-4 are real distributions of dopant for different values of annealing 
time. Increasing of number of curve corresponds to increasing of annealing time 
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Fig.6. Dependences of dimensionless optimal annealing time for doping by diffusion, which have been 
obtained by minimization of mean-squared error, on several parameters. Curve 1 is the dependence of 
dimensionless optimal annealing time on the relation a/L and €= y= 0 for equal to each other values 
of dopant diffusion coefficient in all parts of heterostructure. Curve 2 is the dependence of 
dimensionless optimal annealing time on value of parameter efor a/L=1/2 and €= y= 0. Curve 3 is the 
dependence of dimensionless optimal annealing time on value of parameter € for a/L=1/2 and e€= y= 


0. Curve 4 is the dependence of dimensionless 
a/L=1/2 and e= €=0 
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Fig.7. Dependences of dimensionless optimal annealing time for doping by ion implantation, which 
have been obtained by minimization of mean-squared error, on several parameters. Curve 1 is the 
dependence of dimensionless optimal annealing time on the relation a/L and €= y= 0 for equal to each 
other values of dopant diffusion coefficient in all parts of heterostructure. Curve 2 is the dependence 
of dimensionless optimal annealing time on value of parameter ¢ for a/L=1/2 and €= y= 0. Curve 3 is 
the dependence of dimensionless optimal annealing time on value of parameter & for a/L=1/2 and ¢= 
y= 0. Curve 4 is the dependence of dimensionless optimal annealing time on value of parameter y for 


a/L=1/2 and e= €=0 
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Fig. 8 Normalized dependences of component u, of displacement vector on coordinate z for nonporous 
(curve 1) and porous (curve 2) epitaxial layers 


Further, we have analyzed influence of 
relaxation of mechanical stress on distribution 
of dopant in doped areas of heterostructure. 
Under following condition a< 0 one can find 
compression of distribution of concentration of 
dopant near interface between materials of 
heterostructure. Contrary (at 4>0) one can find 
spreading of distribution of concentration of 
dopant in this area. This changing of 
distribution of concentration of dopant could be 
at least partially compensated by using laser 
annealing [37]. This type of annealing gives us 
possibility to accelerate diffusion of dopant and 
another process in annealed area due to 
inhomogenous distribution, of temperature and 


V(z) 92 


0.0 
0.0 


Arrhenius law. Accounting relaxation of 
mismatch-induced stress in heterostructure 
could leads to changing of optimal values of 
annealing time. At the same time modification of 
porosity gives us possibility to decrease value of 
mechanical stress. On the one hand mismatch- 
induced stress could be used to increase density 
of elements of integrated circuits. On the other 
hand could leads to generation dislocations of 
the discrepancy. Figs. 8 and 9 show 
distributions of concentration of vacancies in 
porous materials and component of 
displacement vector, which is perpendicular to 
interface between layers of heterostructure. 


Zz 


Fig. 9. Normalized dependences of vacancy concentrations on coordinate z in unstressed (curve 1) and 
stressed (curve 2) epitaxial layers 
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CONCLUSION [8] 
In this paper we introduced an approach to 

increase integration rate of field-effect 

heterotransistors in the framework of a 

conventional voltage controlled oscillator. We 

formulate recommendations for optimization of 

annealing to decrease dimensions of the [9] 

considered heterotransistors and to increase 

their density. At the same time one can decrease 

transition area between doped and undoped 

areas of elements of the considered oscillator. 

We formulate recommendations to decrease 

mismatch-induced stress in the considered 

heterostructures. We introduce an analytical [10] 

approach to model diffusion and ion types of 

doping with account concurrent changing of 

parameters in space and time. At the same time 

the approach gives us possibility to take into 

account nonlinearity of considered processes. [11] 
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